PROCEEDINGS 

°F SCIENCE 



00 

o 
o 



DESY 07-104 
SFB/CPP-07-39 
HEPTOOLS 07-017 



(N 
> 

in 

cn 

i> 
o 
l> 
o 



X 



Automatizing tlie application of l\/lellin-Barnes 
representations for Feynman integrals 



J. Gluza^ F. Haas^ K. Kajda" and T. Riemann*^ 

"Department of Field Theory and Particle Physics, Institute of Physics, 

University of Silesia, Uniwersytecka 4, PL-40-007 Katowice, Poland 
^ Deutsches Elektronen-Synchrotron, DESY, Platanenallee 6, 15738 Zeuthen, Germany 



Email: gluza@us.edu.pl Felix.Haas@desy.de kka jda@us . edu .pi 



Tord . Riemann@desy . de 



Feynman diagrams may be evaluated by Mellin-Barnes representations of their Feynman param- 
eter integrals in d = 4- — 2e dimensions. Recently, the Mathematica toolkit AMBRE has been 
developed for the automatic derivation of such representations with a loop-by-loop approach. We 
describe the package and exemplify its use with the e-expansion of the massive one-loop QED 
vertex function. 



XI International Workshop on Advanced Computing and Analysis Techniques in Physics Research, ACAT, 
Amsterdam, The Netherlands, April 23-27, 2007 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



AMBRE - MB-representations for Feynman integrals 



T. Riemann 



1. Introduction 

For many of the precision predictions of observables for LHC and ILC we need the automated 
evaluation of multi-leg and/or multi-loop Feynman diagrams; a summary of the needs and the 
present status may be found in 

A promising approach is the representation of the Feynman diagrams by Feynman parameter 
integrals and the subsequent use of Mellin-Barnes (MB) representations for their further evaluation. 
The method was invented in 1975 for the finite massive scalar three-point function in [^. Massive 
one-loop self-energies and vertices were investigated in d dimensions in An important step 
was the explicit evaluation, in terms of poly logarithms, of the planar on-shell double box in [@-0] 
and of the non-planar case in [^, ^. Slightly different algorithms are used for the derivation of 
proper MB-representations for divergent diagrams as expansions in powers of e = {A — d) /2. The 
algorithm of Tausk was automated in Mathematica and Maple ||^ and in Mathematica [ jlOl ] ; the 
latter package, MB, is publicly available. As input it needs some MB -integral representation and 
performs the e expansion. Quite recently, the Mathematica package AMBRE was published which 



prepares for a large class of Feynman diagrams the MB-integral representation [ 1 1 1. 

In this contribution, we describe the package AMBRE and give a simple example of its use ^ 

2. Construction of Mellin-Barnes representations 

We will use the MB -representation 
1 B 



{A + By 2niT{ 



— - j JaA^B-^r(-a)r(v + a), (2.1) 



where the integration contour separates the poles of the F-functions. 

We evaluate L-loop Feynman integrals^ m d = A — 2e dimensions with N internal lines with 
momenta and masses ni,-, and E external legs with momenta pe'. 

1 [ d''ki...d''kL T{k) 

The numerator T{k) is a tensor in the integration variables: 

T{k) = l,k^,k^kl,... (23) 
The momenta of the denominator functions dj are: 

/ L ^ 

di = qj-mj = = ocnk - PiePe - mj. (2.4) 

\l=l e=l J 

The momentum integrals are replaced by standard Feynman parameter integrals: 

(-l)^-r(jVv-fL) f^^^ ^ ^(^)A^.-.(L+l)/2 



'Examples of different complexity may be found at the webpages | [l^ , |l3[. 

^Often one uses the additional normalization e^^^^; we leave this to the later evaluation with the package MB [flC 
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(2.5) 

with 

N 

A^v = £v,-. (2.6) 

(=1 

The two functions U and F may be derived from 

N 

J/ = Y^Xi{q}-mj) = kMk-2kQ + J, (2.7) 

i=l 

where M/// = l^^j awaaXi, and Qi = Y^Li ocuPiXi, and J = {P?- - m])xu namely: 

U{x) = det(M), (2.8) 
F{x) = -det{M)J + QMQ. (2.9) 

Where we may assume M+ = M. The ?7 and F as well as M = det(M) M^' are polynomials in ;<;. 
We will evaluate L-loop integrals with the loop-by-loop iteration procedure, because the formulae 
simplify for one-loop integrals: 

N 

U = M = M = det(M) = = 1, (2.10) 

i 

N N 

F = -UJ + Q^ = Y,[PiPj-P? + mj]xiXj = Y^fijXiXj. (2.11) 



The simplest tensor factors /"i (r) in (2.5) become: 

= 1, (2.12) 
P,{n = l^xiP^, (2.13) 



p^ik^kP) = i^x^prL-jpf - ^ v- ^^-''^ 



The are the so-called chords introduced in (2^). The general case is: 

Gi(r^) = Giik^'-'-k^-) 

(-l)^v N m r(A^,-£^) 



/j» " '" A; " ' ' 

n^x.xr'-5(i - x^y) I , r^^^^^i {^rP-} 



_rl (Ml. -.Mm 



(2.15) 



with the abbreviations F = F{x) and P = P\{k^) = Y^i^iPi = lLi,e^iPiePe ■ The r starts from zero 
(with i2^^ = 1), and it is =2/^ = for r odd, and j?/^ = g^n^n • ..g^<r-\^ir for r even. The convention 
(/I;, . . .) means the totally symmetric combination of the arguments. 

In AMBRE, the tensorial numerators are assumed to be contracted with a tensor P{m) composed of 
external momenta pe, so that the following quantity is evaluated: 

p(m) Gi(r„,) ^ {p^;---p'^::) G,{k^^---k^-). (2.16) 
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One now has to perform the x-integrations. We do this by the following simple formula: 



1 N 



T- 



i=\ 



r(^i)---r(^jv) 

r(^i H VqN)' 



(2.17) 



From the above text it is evident that the integrand of (2.5) contains besides simple sums of mono- 
mials Y\ix1' (like e.g. Q and M) also sums of monomials with non-integer exponents. This is due 
to the appearance of the factors U{x) and F{x). One may rewrite F{x) and U{x) so that ( |2.17D 
becomes applicable. For the one-loop case this has to be done only for F{x). In (|2.1 1 ), F{x) is 
written as a sum of Np < \N{N + 1) non- vanishing, bihnear terms in xf": 



Nf 



Y^fn{iJ)XiXj 



-(Ny-dL/2) 



(2.18) 



Here, fn{i,j) = fij if fij 7^ 0- Inserting ( |2.18| ) and the tensor function P{T) into ( [2.2| ) allows to 
apply ( 2.17[ ) for an evaluation of the x-integrations. 

As a result, L-loop scalar Feynman integrals may be represented by a single multi-dimensional 
MB -integral and tensor Feynman integrals by finite sums of MB -integrals. 



3. The Mathematica package AMBRE 

In this section we describe the use of the package AMBRE JTl]]. AMBRE stands for Automatic 
Mellin-Barnes Representation. It is a (semi-)automatic procedure written for multi-loop calcula- 
tions. The package works with Mathematica 5.0 and later versions of it. One has to perform the 
following tasks: 

(i) define kinematical invariants which depend on the external momenta; 

(ii) make a decision about the order in which L one-loop subloops (L > 1 ) will be worked out 
sequentially; 

(iii) construct a Feynman integral for the chosen subloop and perform manipulations on the cor- 
responding F-polynomial to make it optimal for later use of the MB representations; 



(iv) use equation ( |2.18| ); 



(v) perform the integrations over Feynman parameters with equation ( [2.17[ ); 



'if useful, one may also consider F(x) with linear and bilinear terms in the Xj. 
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(vi) if needed, go back to step (iii) and repeat the steps for the next subloop until F in the last, 
subloop will be changed into an MB-integral. 

The steps (ii) and (iii) must be analyzed carefully, because there exists some freedom of choice 
on the order of loop integrations in step (ii) and also on the order of MB integrations in step (iii). 
Different choices may lead to different forms of MB-representations. 
The basic functions of AMBRE are: 

• Fullintegral[{numerator},{propagators},{intemal momenta}] - is the basic function for 
input Feynman integrals 

• invariants - is a list of invariants, e.g. invariants = {pl*pl s) 

• IntPart[iteration] - prepares a subintegral for a given internal momentum by collecting the 
related numerator, propagators, integration momentum 

• Subloop[integral] - determines for the selected subintegral the JJ and F polynomials and an 
MB -representation 

• ARint[result,i_] - displays the MB-representation number i for Feynman integrals with nu- 
merators 

• Fauto[0] - allows user specified modifications of the F polynomial fupc 

• BarnesLemma[repr,l,Shifts->True] - function tries to apply Barnes' first lemma to a 
given MB-representation; when Shifts- >True is set, AMBRE will try a simplifying shift 
of variables; the default is Shifts->False 

BamesLemma[repr,2,Shifts->True] - function tries to apply Barnes' second lemma 

4. V312m: the one-loop massive QED vertex 

As an example, we evaluate the one-loop massive QED vertex function V312m (vertex with 3 
internal lines, 2 of them being massive, = 1): 



The corresponding definition in AMBRE: 

» V312m = Fullintegral[{l},{PR[k, 0, nO] PR[k + pi, 1, nl] PR[k - p2, 1, n2]},{k}] 
In a next step we define the invariants: 

» invariants = {pl^- > \,p2^- > \,p\p2- > {s-2)/2} 
With IntPart[l] and SubLoop[integral], we determine the F-polynomial for the diagram. 



As usual, we have to replace here a positive j by j -I- ie. After setting {nl -> 1, n2 -> 1, nO -> 1} 
and applying Barnes' first lemma. 




(4.1) 



F = (xi +X2)^ + [-s]xiX2, 



(4.2) 
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» MBV312m = BarnesLemma[V312m, 1] 

an MB-representation for the Feynman integral is obtained: 



For the subsequent step, the derivation of the e-expansion, the package MB may be used: 

V312m[-l,y] , , , , 

V312m[y] = ^ ^ + V312m[0,y] + £ V312m[l,y] + • • • , 

where we already introduced the conformal variable 



(4.3) 



(4.4) 



(4.5) 



Representations for the first terms of the e-expansion are easily obtained with the following MB 
commands 

» rules = MBoptimizedRuIes[MBV312m, eps -> 0, {}, {eps}] 
» integrals = MBcontinue[MBV312m, eps -> 0, rules] 
» expV312m = MBexpand[integrals, Exp[eps*EulerGamma], {eps, 0, nj] 
We reproduce a few of them here: 



V312m[-l,y] = ;r— / dz{-s) 
V312m[0,y] = ^" dz{-s 



-i-, r3[-z]r[i+z] 
n-2z] 
-1-, r3[-z]r[i+z] 



1 



[Ye - ln{-s) + m>[-2z] - 2^[-z] + »P[1 + z]] ., 
1/4 r+i^+u r3[-z]r[l+z] 



2 

V312m,y] = ^y_^.^^^..(-.) 



(4.6) 



(4.7) 



V312m[2,y] = 



y^+Log[-sf + Log[-s\{-^2YE-^'i'[-2z] +4»F[-z] -2»F[1 +z\) 
+7£(4»P[-2z] -4»P[-z] +2»F[1 +z]) 
-A^>[\ , -2z] + 2»F[1 , -z] + , 1 + z] 
+4(»F[-2z]2 - 2»P[-2z]»P[-z] + + »F[-2z]»P[l + z] 

1/12 r3H]r[i+z] 

2ni J-i^+u ^ ' n-2z] 

a{z) + ^(z)»F(0, 1 + z) + »P(2, 1 + z) + 2»F(0, 1 + z)^ + »I'(0, 1 + z)^ 



(4.8) 



+ 3»i'(0, 1 + z)»I'( 1 , 1 + z) + c(z) [»!'( 1 , 1 + z) + 2»i'(0, 1 + zf] , 



(4.9) 



with some coefficients a(z),&(z),c(z), which depend on Log[—s\,yE,^{k,—z),^{k,—2z) (k = 
0, 1,2). Here, ^[z] and *I'[A:,z] are polygamma functions: 

r'(z) 



PolyGamma [z] = PolyGamma [0,z] = *F(z) 



r(z) 



(4.10) 
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n 2 

PolyGamma [0,n+l] = 5i («) - = ^ - - /£, 

1=1 ' 

PolyGamma [n,z] = ^>[n,z] = = L^^^^- (4-11) 

The integration path includes the straight line ranging from (— /oo + 1 /2) to {+ioo + 1 /2). Closing 
the path to the left, allows to express the integral as an infinite series of residua arising from the 
poles of 

r[\+z],nkA+z], (4.12) 
at z = —n,n = 1,2, •• • with different weight functions G{z), e.g. for V312m [-1 , y ] with 

The resulting respresentations are inverse binomial sums: 

V312m[-l,y] = l£ ^" (4.14) 

V312m[0,y] = - £ ' ^.M ^^' (4.15) 
2„=o ( J(2?i+1) 

V312m[l,y] = \t nn.'",, Hnf + ^2- S2{n)] , (4.16) 
4„=0 ( J(2«+l) 

V312m[2,y] = I (2»^ '^n + 1) [ 1^^^ ~ i'^^ ^"^^^ ^ i ^"^^ ^ ~ ^ 



(4.17) 



etc., and the harmonic numbers Sk{n) are 



" 1 

HarmonicNumber[n,k] = Sk{n) = J^T" (4-18) 

i=i ' 

The simplest of the sums may be done with Mathematica: 

The other sums appearing above may be obtained from sums listed in Table 1 of Appendix D in 
[0: 

°° s" y 

y. 7T1 = -r-^^Hy), (4-20) 



^o Hfi^^+l) ^'^"^ = ^[-4Li2(-y)-41n(y)ln(l+3;) +In2(3;)-2C2] , (4.21) 

y 



Si{n) 



165i,2(-j) - 8Li3(-j) + 16Li2(-3')ln(l +y) 
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-81ii2(l +3;)ln(3;) -41n(l +y)\n^{y) + -ln-'(3;) + 8C2ln(l +y) 



-S2{n) 



-4C2ln(3')-8C3 

y 



toC:){2n+l) 



SUn) 



3(/-l 

y 



(4.22) 
(4.23) 



1 



96Si^2{-y)H^ + j) -965i,3(-j) +4852,2 (-j) 



-24C2ln^(l +j) -481n2(l +3;)Li2(-3;) +48C3ln(l +y) 
+481n(l +y)U3i-y) - 161n(3;)ln^(l +3;) + 24C2ln(3')ln(l +3^) 

+ 121n2(3;)ln2(l +y) -2ln' {y)ln{\ +y) + ^^In'Cj) - 3C2ln'(3') 

+61n2(3;)Li2(-3') +21n2(3;)Li2(3') - 10 l^My) - 24lniy)Usi-y) 



-81n(3;)Li3(y) + 3^4 + 24Li4(-3') + ULUiy) 



/2n 



n=0 



„)(2"+l 



y 



y^-l 



^ln4(3^) + 12Li4(3')+21n2(3;)Li2(3;) 



„foO(2«+l) 



Si{n)S2{n) 



-4C3ln(y)-81n(3;)Li3(3')-12C4 

y 



(4.24) 



(4.25) 



hn\y)\n{\+y)-Un\y)+i;2\nHy) 



+2 In^ (y)Li2 (-y ) + 2 In^ (3;)Li2 (y) + 2^3 ln(3') - 8 \n{y)U^{-y) 



-81n(3;)Li3(3')+2C4 + 16Li4(-y)+12Li4(y) 



(4.26) 



The one-loop QED vertex function is a relatively simple one-scale problem, and in fact one may 
derive the general term of the £-expansion. So far, we applied the Mathematica packages without 
an interference by the user. The packages have been created for the automatization of this type of 
calculations. But here we have an instructive example of the limitations of that. Let us go back to 
Equation ( |4.3[ ) and shift the integration variable z according to z = z' — £, with a related shift of the 
real part of the integration path, u' = u + e. The resulting MB-integral is: 



V312m[y] 



, / d7!{-sr'-' 



r2(-z')r(-z' + e)r(i+z') 



T{l-2e)2niJ ' r(-2z') 
After taking residua as before, but now without using MB, one has to evaluate: 

„er£ ^ r(?i + l + e) 



(4.27) 



V312nifyl = -— V -y- 



(4.28) 



We may apply here the following well-known relation 1 15]^: 
T{n + l+e) 



T{n+\] 



r(l + £)exp 



k=i 



(4.29) 



Express T{—n + e) by r(e). Thereby factors l/(— « + e + fc) are collected, inverse them for small e. 
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and obtain the general term of the inverse binomial sum for the n''^ term in the £-expansion of 



V312m [y] . The first terms agree with ( [^■14[ )-( |4.17| ). The evaluation of the sums may be per- 
formed for the general case by introducing integral representations for the F functions and the 



harmonic numbers (for the latter see e.g. [|16|]), and then performing the two-dimensional integral. 

Alternatively to the approach described here, one may determine the vertex function in terms 
of Harmonic Polylogarithms (HPLs) JXt] ] by the differential equation method, which allows easily 
to derive the £ -expansion of the Feynman integral SE2 12m. The vertex may be related to SE2 12m 
and to the tadpole Tlllmby integration by parts; see e.g. [Ill 



V312m : 



{d - 2)Tlllm - l{d - 3)SE212m 
(^-4)(^-4) ■ 



(4.30) 



The master integrals Tlllm and SE212m are tabulated e.g. in the file masters HP L.m [ ]19 , 
20]. From the corresponding formulae it may be seen that only HPLs with arguments and -1 
appear. An approach to describe the vertex in d dimensions with hypergeometric functions and the 



subsequent derivation of the £-expansion may be found in |21]. 



5. Summary 

To summarize, for many applications of present phenomenological or more theoretical interest 
the package AMBRE solves an important part of the calculational problem: the semi-automatic 
derivation of MB-representations for a large class of Feynman integrals. 
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